Critical behavior of the 3d random field Ising model: 
Two-exponent scaling or first order phase transition? 
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In extensive Monte Carlo simulations the phase transition of the random field Ising model in 
three dimensions is investigated. The values of the critical exponents are determined via finite size 
scaling. For a Gaussian distribution of the random fields it is found that the correlation length £ 
diverges with an exponent v = 1.1 ± 0.2 at the critical temperature and that x ~ £ 2- '' with r\ — 
0.50±0.05 for the connected susceptibility and Xdis ~ ^ 4 ^ r > with rj = 1.03±0.05 for the disconnected 
susceptibility. Together with the amplitude ratio A = lmiT^T c Xdis/x 2 (hr/T) 2 being close to one 
this gives further support for a two exponent scaling scenario implying rj = 2n. The magnetization 
behaves discontinuously at the transition, i.e. = 0, indicating a first order transition. However, 
no divergence for the specific heat and in particular no latent heat is found. Also the probability 
distribution of the magnetization does not show a multi-peak structure that is characteristic for the 
phase-coexistence at first order phase transition points. 

75.10.H, 05.50, 64.60.C 



^1- 

o 
m 
O 
in 
o\ 
+-> 



c 

o 
o 



X 



I. INTRODUCTION 



The critical behavior of the random field Ising model 
(RFIM) has been investigated intensively over two 
decades now (seed for a review) but a full agreement 
over the most fundamental issues still missing. Although 
it has been shown rigorouslyH that at low temperatures 
and sufficiently small field strength there is indeed long 
range ferromagnetic order in three dimensions, the char- 
acterization of the phase transition is still unclear. The 
idea of dimensional reductionEl failed to describe the crit- 
ical behavior of the model correctly. The reason for this 
was essentially the presence of many local minima in the 
free energy landscape that do not allow for a straight- 
forward perturbative treatment. A reminiscence of di- 
mensional reductioruis still present in the scaling theo- 
ries for the RFIMta where the dimension d has been 
replaced by d — 9 in hyper-scaling relations and 9 was 
introduced as a third exponent independent of v and 
r\. However, these (droplet)-theories relied upon the as- 
sumption of a second order phase transition. Part of the 
results obtained via Monte Carlo simulations support this 
scenaricoJ, whereas other simulations find indications for 
the transition to be of first ordeialld. Moxe surprisingly 
no latent heat was found at the transition!!!], contrary-to 
what one might expect at a first order phase transition!!!!, 
where a discontinuity in the internal energy at the critical 
point causes the specific heat to diverge with the volume 
of the system. Actually no .-divergence at all could be 
detected for the specific heatE3. 

The experimental situation is also far from clear. Ex- 
periments on diluted antiferromagnets in an uniform ex- 
ternal field (DAFF), which have been demonstrated to 
be in the same universality class as the RFIMI3L3, pro- 



vided evidence for a second order phase transition with 
varying jCstimai.es for the critical exponents (for a discus- 
sion seeEl). IrJiil a logarithmic singularity for the specific 
heat was found, implying a = XL Quite recently, subse- 
quent to the prediction made inE3, experimental evidence 
of a cusp-like-singularity of the specific heat, i.e. a < 0, 
was reportedE-3. Moreover, other possible scenarios, like 
a first order traasitiontj or a so called "trompe l'oeil crit- 
ical behavior"l!H have also been discussed. Experiments 
addressing the equilibrium quantities of the RFIM (or the 
DAFF, respectively) are extremely difficult to perform, 
since by approaching the transition the system falls out 
of equilibrium veoz-rapidly and freezes into a metastable 
domain state (seeE3 for a review) . Therefore these issues 
have not been settled by now. 

The purpose of this paper is to estimate all critical ex- 
ponents with the help of extensive Monte Carlo investiga- 
tions of the RFIM in three dimensions. In a recent paper 
results for the exponent&jaf the RFIM with a binary dis- 
tribution were presented!!!! In mean field theory!!!! there 
is a sharp discrepancy between the phase diagram for a 
binary distribution and that of a continuous distribution 
with a finite probability for zero random fields: in the 
former case the phase diagram has a tritcritical point, 
where the phase transition changes from second order 
for low fields to first order at higher fields; in the latter 
case the transition is always of second order. One might 
speculate that a similar scenario occurs also-Ln three di- 
mensions and that the indications found in£3 for a first 
order transition are in fact an artifact of the binary dis- 
tribution. Therefore it is important to perform a similar 
analysis for a continuous distribution, which is presented 
in this paper. 

It is organized as follows: In section 2 we present the 
model and review various scaling predictions that are 
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supposed to hold in case he transition is of second or- 
der. Section 3 reports our results of the finite size scaling 
analysis and section 4 discusses the possibility of a first 
order transition. Section 5 is a summarize and an out- 
look. 



II. MODEL AND SCALING THEORY 

The random field fsing model (RFIM) in three dimen- 
sions is defined by the Hamiltonian 



H 



(i) 



(id) 



where Si = ±1 are Ising spins, the first sum is over near- 
est neighbor pairs on a simple cubic lattice and the hi 
are independent quenched random variables with mean 
zero [hi] av — and variance [hf] av = h%- ft has been 
shown rigorouslyo that this model has a transition to 
long range ferromagnetic order along a line h r (T) in the 
(h r , T)-diagram. 

If this transition is second order, a scaling theory has 
been proposecuQ that relies upon the assumption that 
random field induced fluctuations dominate over thermal 
fluctuations at T c . This implies for the singular part of 
the free energy 



F ■ 

± sing 



(2) 



where £ is the correlation length diverging at the critical 
temperature like 



with t = T — T r 



(3) 



and 9 is a new exponent. Random field fluctuations 
alone produce typically an excess field of the order of £ d / 2 
within a correlation volume, so a naive guess would be 9 
roughly 1.5 in three dimensions. From (^|) the modified 
hyper-scaling relation follows 



2 — a — v{d 



(4) 



a being the specific heat exponent. It also implies an 
exponential divergence of the relaxation time r at T c : 
t ~ exp(A£ e ). The decay of the connected and discon- 
nected correlation functions at T r 



C(v) = [{S S r )-(S )(S r )] 3 
Cdis(r)= [<So)<S r )] av 



-{d-2+v) 
-(d-4+5f) 



defines two exponents 77 and 77. These are expecte 
be related to the new exponent 9 via 

9 = 2 - 77 + 77 . 



(5) 
(6) 

to 
(7) 



By approaching the critical temperature from above the 
corresponding connected and disconnected susceptibili- 
ties diverge as 



X = JdrC(r) ~ (T-T c )-< 
Xdis= /drC dis (r)~ (T-T c )-~ 

These exponents are related by 

7 = i/(2 — r]) , 
7 = i/(4 - 77) . 



(8) 
(9) 



(10) 
(11) 



In addition one has the usual scaling relations for the 
magnetization exponent 



(3=(d-4 + rj), 
and the Rusbrooke equality 

a + 2/3 + 7 = 2 . 



(12) 



(13) 



Obviously there seem to be three independent critical ex- 
ponents and a central issue of the activities on the critical 
properties of the RFIM is the quest fcpan additional scal- 
ing relation. Already Imry and Mao conjectured that 
.Fsing ~ X, where \ ~ (T - T c )~ 
This would imply 9 = 2 — 77 via 



7 is the susceptibility. 
(0) and therefore with 



(14) 



77 = 2? ? 



A set of analytical arguments by Schwartz et alBl 
supports this two-exponent scaling scenario indicated, 
by (|l4|). Hence the exact Schwartz-Soffer inequalitytil 
77 < 277 might be fulfilled as an equality, as in,-ffl4). A 
recent high temperature series analysis providedostrong 
evidence of this equality to hold by estimating 7 and 7 
via Pade-approximation. 

In this paper we intend to determine all critical ex- 
ponents with Monte Carlo simulations. The model is 
defined by the Hamiltonian on a finite lattice of lin- 
ear size L with periodic boundarv-conditions. Instead of 
a binary distribution, as used inE£!l, we consider here a 
Gaussian distribution 



P{hi) 



eM-h 2 j2hl) 



(15) 



for the random fields. Unfortupately this makes the use 
of multispin coding impossiblccj, which slows down the 
spin-update algorithm (Metropolis) by approximately a 
factor 10. The simulations were done at a constant ratio 
h r /T = 0.35. 

As was pointed out iriiJ it is of paramount importance 
to perform an extensive disorder average, since quanti- 
ties like the susceptibility \ are non-self-averaging in this 
model. Therefore at least 1000 samples for each tem- 
perature and system size were used. Equilibration was 
checked by starting each sample with two different ini- 
tial configurations, one with all spins down and the other 
with all spin up, and running as long as both "replica" 
yield the same values for all thermodynamic quantities 
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of interest. For L = 16 up to 2 ■ 10 6 Monte Carlo sweeps 
were necessary to fulfill this criterion. 

The following quantities were calculated for each disor- 
der realization: The average magnetization per spin (M), 
its square (M 2 ), the average energy per spin (E) and its 
square (E 2 ). The angular brackets, (. . .), denote a ther- 
mal average for a single random field configuration and 
M = E = H (S)/N and N = L d . From 

these data we get the specific heat per spin, C, the sus- 
ceptibility Xi the- disconnected susceptibility, Xdis, and 
the order parameter, m, as follows 

[C] av = N{[(E 2 )} av - [<£} 2 ] av }/T 2 , (16) 

[m} av = [\(M)\] av , (17) 

[ x ] av = 7v{[(Af 2 )] av -[(Af) 2 ] av }/T, (18) 

[Xdis]av = N[(M) 2 ] av , (19) 

where [. . .] av denotes the average over different random 
field configurations. The finite size scaling function of 
any quantity A(T, L), which behaves in the infinite sys- 
tem like Aqo ~ t~ a by approaching the critical tempera- 
ture, reads 



A(T,L) = L a / v A{tL 1 / v ) 



(20) 



In case of the specific heat the scaling function C(x) has a 
maximum at some value x = x* . For each lattice-size we 
estimate the temperature T*(L), where T 2 [C] av is max- 
imal. Since t*{V)L x l v = x* we obtain in this way the 
critical temperature T c and the correlation length expo- 
nent v from 



t*(L) = T*(L)-T c = x*L 



-l/v 



(21) 



We denote the value of T*(L) 2 [C] av at this temperature 
T*(L) by [C*1 av and similarly for the other quantities in 



Eq. (16 



T*(L)[ X (L,T*(L))U, [ X * i( 



[ X dis(L,T%L))] av and [m*] av = [m{L,T*{L)% 
one expects the following behavior: 



[C 

Ix 



I av 



lsjav 



L a ' v C{x*) , 
L 2 ^_ X (x*) , 

L^xaUx*) , 
L-Pl v m(x*) . 



Then 

(22) 
(23) 
(24) 
(25) 



Through the size dependence of [C*] av , [x*]av, [Xd is ]av 
and [m*] av we obtain the exponents a, rj, rj and (3. In 
the vicinity of x* the scaling function C(x) can be ap- 
proximated by a parabola. Therefore three temperatures 
near the maximum of the specific heat are enough to de- 
termine the values of T*{L) as well as [C*] av etc. Let us 
assume that we measured C% ± Ai, C 2 ± A 2 and C3 ± A 3 
for temperatures T\ < T2 < X3. In order to get a reli- 
able estimate for C(L) as well as its errorbar we have to 
have Ci + Ai < C 2 - A 2 and C 3 + A 3 < C 2 - A 2 , hence 
it is important to have very small statistical errors Ai, 



A 2 and A3. The difference between the temperatures 
should be small, otherwise the second order polynomial 
approximation fails, but not too small, since then the 
above requirement for the difference between C\ , C 2 and 
C 3 cannot be fulfilled. We used T 2 - T x = T 3 - T 2 = 0.1 
for the small lattice sizes (L < 10) and 0.05 for the larger 
sizes. Once T*(L) is determined, the value of the quan- 
tities in (p2]]25| ) are again determined by fitting a second 
order polynomial. Finally, the exponents are determined 
via least square fits of the data obtained to the functional 
forms given in (I22H2C 



III. THE CRITICAL EXPONENTS 

First we determine the temperature T*(L), where the 
maximum of the specific heat occurs and fit the corre- 
lation length exponent v in such a way that the data 
for T*(L) lie on a straight line if plotted against L^ 1 ^. 
Then relation ([H]) is fulfilled and T c can be read of from 
the intersection of this line with the y-axis. The result is 
depicted in figure (|l|), one obtains 



= 1.1 ±0.2 and T r = 3.695 ± 0.02 



(26) 



Next we analyse the size-dependence of maximum of 
the specific heat. According to ( p2| ) one should be able 
to read off the exponent-ratio a/v from the slope of a 
straight line fit in a log- log plot of [C*] av versus L. The 
first observation is that this procedure does not work 
since the data points show a significant negative curva- 
ture in a log- log plot, indicating a much weaker depen- 
dency then algebraic. In reminiscence of the experimen- 
tal finding of a logarithmic divergence of the specific heat 
(i.e. a — 0) we plotted [C*] av versus log(L), and still 
found a (now less pronounced) negative curvature in the 
data points, indicating a still weaker or even no diver- 
gence. Finally we hypothesized a cusp-like singularity, 



(27) 



with a < 0. The result of the corresponding fitting pro- 
cedure is shown in figure [| As one can see the data are 
fully consistent with this type of dependency and we get 



u/i> 



-0.45 ±0.05 and 



30.4 ±0.1. 



(28) 



With the estimate for v from (|26|) it is a = —0.5 ± 0.2. 
Although the ratio a/v can be determined quite accu- 
rately, the value for a alone depends on v which has a 
much larger large error bar. 

According to (^) the susceptibility at the temperature 
T*(L) should diverge algebraically, and the exponent r] 
is obtained by a least square straight line fit of [x*]av 
versus L in a log-log plot, which is depicted in figure ||. 
The result is 



j] = 0.50 ±0.05 



(29) 
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from which one obtains via jic| ) 7 = 1.7 ± 0.2. 

The finite size scaling behavior of x near the critical 
temperature is expected to be 



[ X ] av (T,L)«L 2 -"x(^(T-T c )) 



(30) 



In figure ^ the corresponding finite size scaling plot using 
the values for T c , v and r\ from (|2^) and ( p9|) is shown. 

The disconnected susceptibility at T* (L) is depicted in 
figure ||, which should diverge according to ( p3 ) with sys- 
tem size. Via a least square fit of the data to a straight 
line in the log-log plot one obtains the estimate 

rj= 1.03 ±0.05 . (31) 

With the scaling relation (0) one can estimate 6 to be 

9 = 1.5 ±0.1, (32) 

which is in good agreement with the value 8 — d/2 men- 
tioned in the last section. 

Comparing our values for rj ( |2S| ) and rj (^lj) one notes 
that the relation ([l4|) might be fulfilled, which would sup- 
port the two-exponent scaling scenario mentioned in the 
last section. In order to check this we calculated the 
amplitude ratio 



A = lim 



[Xdh 



(33) 



as suggested where in an extensive high- 

temperature analysis it was found that A — 1, implying 
strong evidence for two-exponent scaling. The finite size 
scaling form of © is A(L, T) = A(L 1 / iy (T - T c )) and 
therefore we expect 



A*{L) = A(L,T*(L)) = 1 



(34) 



independent of system size. The result that we obtain 
from our MC-data is table 1: 



L 


A* 


AA* 


4 


1.19 


0.22 


6 


1.14 


0.23 


8 


1.10 


0.33 


10 


1.12 


0.29 


12 


1.00 


0.26 


16 


0.80 


0.27 



TABLE I: The amplitude ratio A*(L) = Ail., T*{L)\ and its 
estimated error for different system sizes. 

As one can see A* is close to one, which supports the 
equality (14). Because of the large errorbars the slight 
decrease in the estimate of A* with increasing system size 
seems not to be significant. 

From equation ( |3l| ) one concludes that the discon- 
nected susceptibility diverges roughly with the volume 
of the system, since from (11) one obtains 7/1/ = 4 — rj k, 
3 = d. If one looks at equation (|l9|), where Xdis is defined, 



this implies that [(M) 2 ] av = Xdis/L d is then independent 
of system size at the critical temperature. Since the mag- 
netization is defined in a similar manner ( |l7| ) , namely the 
square replaced by the modulus, it does not come as a 
surprise that one observes that also the magnetization is 
also independent of system size at the critical point. This 
connection manifests itself in the relation ( |L2] ) , which pre- 
dicts (3 ~ if rj ~ 1. By looking at [rn*] a v as a function of 
L we cannot detect any significant variation with system 
size. As an even stronger evidence we show in figure [| a 
finite size scaling plot of the magnetization according to 



L 1/v (T~T c )) if /3 = 



(35) 



using the values for T c and v estimated above. Also by 
studying the probability distribution of the magnetiza- 
tion in the next section, we see that that the system in the 
thermodynamic limit is already ordered at T c . Einally, if 
the transition is first order, one expects 77 = l/2c3, which 
is in agreement with our estimate (p9|). 

Finally we study the dimcnsionlcss coupling constant 
or Binder cumulant 



g(T,L) = 0.5^3- 



[(M 2 ) 2 } 



(36) 



which is expected to be independent of size at T c . In fig- 
ure [7] we show the result for three different system sizes 
as a function of temperature. The intersection of the 
curves of the two largest system sizes lie very close to the 
estimated value of T c = 3.695. 

One observes that close to the critical temperature the 
value of the cumulant g lies already very close to 1 . This 
means that the probability distribution for the order pa- 
rameter is centered around a non- vanishing value for the 
magnetization, implying a large degree of magnetic or- 
der already at T c . This is fully compatible with (3 = 
and a discontinuous jump in the expectation value for 
the disorder-averaged magnetization at the critical tem- 
perature 



IV. PROBABILITY DISTRIBUTIONS 

As mentioned above the susceptibility is a highly non- 
selfaveraging quantity, which can be demonstrated by 
looking at the distribution P(x) f° r the probability to 
find a value x f° r f ne susceptibility of one sample (i.e. 
realization of the disorder). P{x) is f ar from being a 
Gaussian distribution and it possesses pronounced long 
tails as shown in figure |H (see alscEa). 

This is one of the predictions of the droplet theorj 
With a probability proportional to L~ e a sample has a 
vanishing thermodynamic expectation value for the mag- 
netization leading to a strongly enhanced susceptibility 
that is in this case equal to the disconnected susceptibil- 
ity, as can be seen from its definition (|l^,|l^) . This argu- 
ment leads to the equation (0). Moreover it says that the 
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second moment [x 2 ]av of the distribution P(x) is domi- 
nated by these rare events occurring with probability L~ e 
and giving a contribution proportional to the square of 
the disconnected susceptibility, i.e. (L 4 ~'') 2 . This gives 

[X 2 U(L,T*(L)) , (37) 

with C = -0+2(4-77) w 4.5 via @ and @. This gives 
a much larger contribution than the square of the mean 
value [x]av ~ L 2{ - 2 ~' qS) « L 3 . By explicitly calculating 
[x 2 ]av(T'*L, L) with our MC-data we get an estimate of 
£ = 4.0 ± 0.2, which is indeed closer to the above droplet 
prediction of 4.5. 

The scaling hypothesis concerning the probability dis- 
tribution of the magnetization can be checked by inspect- 
ing the latter directly. In figure || we show the probabil- 
ity distribution of the magnetization for different system 
sizes at a temperatures close to T* (L) . Let us focus at the 
moment on the behavior of Pi(m) at zero magnetization. 
Note that m is defined as the modulus of the thermody- 
namic expectation value of the fluctuating quantity M 
and can thus take on only positive values. 

At a usual second order phase transition Pl (m, T*(L)) 
would have a peak at m = with a width that is pro- 
portional x*(V)jL d ~ L 2 ~ v ~ d . Thus this peak narrows 
and gains more weight with increasing system size. Here 
it is quite contrary, the peak develops at a nonzero value 
TOqo for the magnetization. If it would be a conventional 
Gaussian distribution centered around this nonvanishing 
most probable value for the magnetization, its value at 
m = would decrease exponentially with system size: 
P *(L) = P L (m = 0,T*(L)) oc exp(-m 2 0O L d - 2+r i) . How- 
ever, according to the above mentioned scaling theory, 
one expects 

P5{L)~L-°, (38) 

which is a much slower decay. In order to improve 
the statistics we estimated Pq(L) via P(m = 0) <w 
Am -1 J Q Am dm P(m) with Am = 0.2 - 0.3. In figure 
[Toj we depict the result from which we conclude that ( |38| ) 
obeyed and we estimate 6 = 1.0 ± 0.1, which is smaller 
than the value (|32|) obtained by the relation It is 
clear that a nonvanishing Am overestimates Pq(L) more 
for larger system sizes, which might explain this discrep- 
ancy. 

Finally we turn our attention to the..va,riation of P(m) 
with temperature. According to ref.c3E3 the probability 
distribution of the order parameter provides a mean to 
discriminate between a first and a second order phase 
transition. At a usual first order transition one expects a 
phase-coexistence at the critical point, which means that 
P(m) has two significant peaks, one at zero (represent- 
ing the high temperature phase) and one at a positive 
value moo > (representing the low temperature phase) . 
In this case the effective transition temperature T C (L) of 
the finite system can be identified b« the temperature, 
where both peaks have equal weightEZl. 



If we look at figure O, where Pt,=v&{jn] is shown for 
various temperature values, we detect no sign of a char- 
acteristic double-peak structure mentioned above. On 
the contrary, there is only a single peak moving continu- 
ously to the right with decreasing temperature. Thus we 
conclude that phase-coexistence at T c can be excluded 
in the RFIM, although the magnetization jumps discon- 
tinuously there. Hence there is a discontinuity in one of 
the derivatives of the free energy (namely with respect 
to a homogeneous external magnetic field), which means 
that the transition is first order in a strict sense. Never- 
theless it is a quite unusual first order phase transition: 
neither a latent heat nor a phase-coexistence at T c could 
be detected in our simulations. 

We would like to remark that our findings for the prob- 
ability distribution Pl (m, T) are fully compatible with, 
the scenario, which is suggested by the droplet theorytrEl 
— in particular the systematic decrease of the probabil- 
ity for samples with zero magnetization J3S|) . Once we 
accept that Pl(tti = 0) shrinks to zero at T c in the limit 
L — > oo the weight of the distribution has to cumulated 
at a nonzero most probable value for m. As we have dis- 
cussed this implies a discontinuous jump of the average 
magnetization at T c . Thus one might speculate that our 
finding j3 = is an intrinsic feature of the droplet theory, 
although such a possibility has not been discussed yet to 
our knowledge. 

Another paradigm of usual first order phase transitions 
is that the correlation length stays finite at T c , although 
it might become very large (for a counter-example seec3). 
In this case the scaling analysis performed in the last sec- 
tion is unjustified, since it relies on a second order phase 
transition scenario with a diverging correlation length. 
However, since neither the specific heat nor the order 
parameter probability distribution behaves as usual at a 
first order transition there seems also to be no reason us 
for to expect a finite correlation length at T c . Finally let 
us mention that we also looked at the probability distri- 
bution of the internal energy. Here no remarkable behav- 
ior could be observed. 



V. DISCUSSION 

To summarize we list the values for the critical expo- 
nents that we have calculated in table II together with 
those for the binary distribution obtained inE2l. Let us 
compare these numbers with other known values from 
the literature: 77 = 0.25 ± 0.03, rj - 0.8, 7 = 1.7 ± 0.2 
and indications of a first order transition from MC- 
simulations of the RFIMH; r) = 0.5i43.1 and rj = 1.0±0.3 
from MC-simulations of the DAFH3; rj = 1.1 ± 0.1 and 
v = 1.0 ± 0.1 from a ground-state investigation via com- 
binatorial optimizationB; 7 = 1.9 —, J2..2 from real space 
renormalization group calculationsc3; 7 = 2.1 ± CU2, 
7 = 4.2 ± 0.4 from high temperature series expansionEj; 
v = 1.25 ± 0.11, rj = 0.89 ± 0.10 and 77 = 0.4 ± 0.8 from 
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a weighted mean field theoryeO. Our estimates are well 
compatible with all of these results. 





Gaussian 




Binary 




'1 


0.50 


±0.05 


0.56 


±0.03 


V 


1.03 


±0.05 


1.00 


±0.06 


e 


1.53 


±0.1 


1.56 


±0.1 


V 


1.1 


±0.2 


1.6 


±0.3 


7 


1.7 


±0.2 


2.3 


±0.3 


7 


3.3 


±0.6 


4.8 


±0.9 


P 


0.00 


±0.05 


0.00 


±0.05 


a 


-0.5 


±0.2 


-1.0 


±0.3 



TABLE II: The critical exponents for the Gaussian distri- 
bution considered in this paper and for the binary distri- 
bution investigated in ref.u, both for the constant field- 
strength/temperature ratio h r /T = 0.35 

The first important observation is that the results ob- 
tained for a Gaussian and a binary distribution do not 
differ significantly. In particular the indications for a 
discontinuity in the magnetization [3 = detected for 
the binary distribution in ref.ta are also present for the 
Gaussian distribution investigated in this paper. This 
is a crucial point with regards to mean-field theoryc3, 
where an essential difference between the two kinds of 
distributions is predicted: In contrast to the continuous 
case with nonvanishing weight at zero field strength the 
binary distribution is expected to have a tricritical point, 
where the transition changes from second to first order. 
Our conclusion in this paper is that in three dimensions 
both distributions lead to identical results for the inves- 
tigated field strength. 

Our estimate for the correlation length exponent v has 
a rather large errorbar, which is also the main source for 
the error in the values for 7, 7 and a. One notes also 
a discrepancy between the values of v for the Gaussian 
and for the binary distribution, but they remain compat- 
ible within the errorbars. The Schwartz-Soffer inequality 
rj < 2rj seems to be fulfilled as an equality, which sup- 
poets the two-exponent scaling scenario, as also found 
ir£3. An explicit calculation of the amplitude ration in 
section III gives further evidence in this respect. 

The analysis of the probability distribution of the sus- 
ceptibility showed pronounced long tails as predicted by 
the scaling theoriestm. The source for these long tails 
are rare samples possessing a vanishing net-field and thus 
having a vanishing expectation value for the magnetiza- 
tion. In fact these rare samples provide the dominant 
contribution to the higher moments of the distribution, 
which we checked explicitly. The analysis of the prob- 
ability distribution of the magnetization supports this 
picture and shows that the probability of these rare sam- 
ples indeed decays like L~ e with system size, where 9 is 
an exponent whose value turns out to be compatible with 
other predictions of the scaling theory. 

It turns out that the magnetization exponent (3 is 
zero, which means that the magnetization jumps dis- 



continuously from zero to a finite value at the criti- 
cal temperature. Although this might be called a first 
order phase transition a closer inspection of the prob- 
ability distribution of the magnetization revealed the 
absence of a double-peak structure characterizing the 
phase-coexistence usually present at a first order phase 
transition. 

Furthermore, at a first order phase transition the 
specific_heat usually diverges with the volume of the 
systemtj, but the exponent a is negative, which means 
that the specific heat does not diverge at 3i™ By means 
of birefringence techniques Belanger et al.E-3 concluded 
from their experiments on Feo.47Zno.53F2 that a = 
(i.e. a logarithmic divec&ence of the specific heat). Only 
recently it was shownlij that the same kind of experi- 
ments on Feo.ssMgo.isB^ are better compatible with a 
cusp like singularity of the specific heat and a = — f , 
concurring with the value reported in table II. However, 
this value for a together with the other estimates in ta- 
ble II would violate the modified hyper-scaling relation 
(Q). Moreover, Schwartzes derived an exact inequality 
2 — a < v{d — 2 + r\) and accepting 7/ w 0.5 (since this 
result has a much smaller errorbar than v) it would im- 
ply v > 2, which is a gather large value. Indeed in a re- 
cent Migdal-KadanofEl such a large exponent was found: 
v = 2.25, and also report a = —1.37 and (3 — 0.02, con- 
sistent with table II. However, a value for a that is nega- 
tive and large in modulus causes serious difficulties with 
respect to the Rushbrooke relation ([l3]) and ±he more 
rigorous Rushbrooke inequality a + 2(3 + 7 > 2c3. 

Let us conclude with the remark that this investigation 
on one hand confirms-the picture of the phase transition 
scenario found in ref.E3. On the other hand it poses some 
serious puzzles concerning various scaling relations for 
the exponents that still remain to be solved in the future. 
In addition quite recently the possibility of an intermedi- 
ate spin glass phase, located between the para magn etic 
and the ferromagnetic phase has been discussedE3E3. It 
seems to us to be quite important to clarify the conse- 
quences of the existence such a phase, which is supposed 
to posses a finite correlation length, for the nature of the 
transition considered here. 



ACKNOWLEDGEMENT 

I would like to thank A. P. Young for his encourage- 
ment and many extremely helpful discussions. The com- 
putations were done on the Intel Paragon System from 
the Supercomputer center (HLRZ) at the Forschungszen- 
trum Julich. This work work was performed within the 
SFB 341 Koln-Aachen- Julich, supported by the DFG. 



6 



1 D. P. Belanger and A. P. Young, J. Magn. Magn. Mat. 100, 
272 (1991), H. Rieger, Monte Carlo Studies of Spin Glasses 
and Random Field Systems in: Ann. Rev. Comp. Phys. II, 
World Scientific, Singapore (1995). 

2 J. Bricmont and A. Kupiainen, Phys. Rev. Lett. 59 1829 
(1987). 

3 A. Aharony, Y. Imry and S. K. Ma, Phys. Rev. Lett. 37, 
1364 (1976); A. P. Young, J. Phys. C 10, L257 (1977); G. 
Parisi and N. Sourlas, Phys. Rev. Lett. 43, 744 (1979). 

4 J. Villain, J. Physique 46, 1843 (1985). 

5 D. S. Fisher, Phys. Rev. Lett. 56, 416 (1986). 

6 A. J. Bray and M. A. Moore, J. Phys. C 18, L927 (1985). 

7 A. T. Ogielski and D. A. Huse, Phys. Rev. Lett. 56, 1298 
(1986). 

8 A. T. Ogielski, Phys. Rev. Lett. 57, 1251 (1986). 

9 A. P. Young and M. Nauenberg, Phys. Rev. Lett. 54, 2429 
(1985). 

10 H. Rieger and A. P. Young, J. Phys. A 26, 5279 (1993). 

11 M. Nauenberg and B. Nienhuis, Phys. Rev. Lett. 33, 944 
(1974). 

12 S. Fishman and A. Aharony, J. Phys. C 12, L729 (1979). 

13 J. Cardy, Phys. Rev. B 29, 505 (1985). 

14 D. P. Belanger, A. R. King and V. Jaccarino, Phys. Rev. 
B 31, 4538 (1985). P. Pollack, W. Kleemann and D. P. 
Belanger, Phys. Rev. B 38, 4773 (1988). 

15 M. Karszewiski, J. Kushauer, Ch. Binek, W. Kleemann and 
D. Bertrand, J. Phys. C 6, L75 (1994). 



R. J. Birgeneau, R. A. Cowley, G. Shirane and H. 

Yoshizawa, Phys. Rev. Lett. 54, 2174 (1985). 

J. P. Hill, Q. Feng, R. J. Birgeneau and T. R. Thurston, 

Phys. Rev. Lett. 70, 3655 (1993); Z. Phys. B 92, 285 

(1993). 

W. Kleemann, Int. J. Mod. Phys. B 7, 2469 (1993). 

A. Aharony, Phys. Rev. B 18, 3318 (1978). 

Y. Imry and S. K. Ma, Phys. Rev. Lett. 35, 1399 (1975). 

M. Schwartz and A. Sofier, Phys. Rev. Lett. 55, 2499 

(1985). 

M. Schwartz, J. Phys. C 18, 135 (1985); M. Schwartz and 
A. Soffer, Phys. Rev. B 33, 2059 (1986); M. Schwartz, M. 
Gofman and T. Nattermann, Physica A 178, 6 (1991). M. 
Schwartz, Europhys. Lett. 15, 777 (1994). 
M. Gofman, J. Adler, A. Aharony, A. B. Harris et al., Phys. 
Rev. Lett. 71, 1569 (1993). 

H. Rieger, J. Stat. Phys. 70, 1063 (1993). 

I. Dayan, M. Schwartz and A. P. Young, J. Phys. A. 26, 
3093 (1993). 

K. Binder and D. P. Landau, Phys. Rev. B 30, 1477 (1984). 
M. S. S. Challa, D. P. Landau and K. Binder, Phys. Rev. 
B 34, 1841 (1986). 

P. W. Anderson and G. Yuval, J. Phys. C 4, 607 (1971). 
M. S. Cao and J. Machta, Phys. Rev. B 48, 3177 (1993). 
D. Lancaster, E. Marinari and G. Parisi, Weig hted Mean 



Field Theory for the RFIM, |cond-mat/9412069 . 
G. S. Rushbrooke, J. Chem. Phys. 39, 842 (1963). 
M. Mezard and A. P. Young, Europhys. Lett. 18, 
(1992). 

M. Mezard and R. 
(1994). 



653 



Monasson, Phys. Rev. B 50, 7199 



7 



0.1 0.2 0.3 

L -i/v 

FIG. 1. The temperature T*(L), where the specific heat attains its maximum, versus L 1//n ", with v as in (pt^). From right 
to left one has L — 4, 5, 6, 8, 10, 12 and 16. The full line is a least square straight line fit. 
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FIG. 3. The value of the susceptibility [x]av at the temperature T*(L) versus system size L in a log-log plot. The full line 
is a least square straight line fit, which gives a slope of 1.55 ± 0.05 that is an estimate for 2 — r/ according to 
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FIG. 4. A scaling plot of the susceptibility [x]av with the parameters for T c , v and 77 given in (^lj) and 
connects only the points for L — 6 as a guide for the eyes. 
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FIG. 5. The disconnected susceptibility [xdis]av at the temperature T*(L) versus L in a log- log plot. The full line is a least 
square straight line fit, which gives a slope of 2.97 ± 0.08 that is an estimate for 4 — 7j according to 
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FIG. 6. A scaling plot of the magnetization [M] av with the parameters for T c and v given in The full line connects 

only the points for L = 6 as a guide for the eyes. Note that the data for [M] av are not rescaled by a factor L~^^ ' , implying 
that indeed (3 — 0. 
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FIG. 7. The cumulant g as defined in (j36|) versus temperature for three different system sizes. The full curves are least 
square fits of third order polynomials to the data points. The intersection of the L — 16 and L = 12 curves lies already close 
to the estimated critical temperature T c = 3.695 (Elh . 
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FIG. 8. The probability distribution P(x) for L = 16 at the temperature T = 3.75 « T*(L = 16). Note the logarithmic 
scale of the abscissa. 
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FIG. 9. The probability distribution Pi(m) for increasing system sizes close to the temperature T*(L). The tendency 
towards a dominant peak at nonvanishing magnetization for increasing system size is obvious. 
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FIG. 10. The probability Pz,(m = 0, T = T*(L)) versus system size L in a log-log plot. The full curve is a least square 
straight line fit with slope 6 — — 1.0±, 0.1. 
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FIG. 11. The probability distribution P(m) for L = 16 at different temperatures. Note that T = 3.80 is slightly above 
T*(L = 16) and T = 3.75 is slighly below it, but still above the estimated value for T c = 3.695. 
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